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ABSTRACT 

We discuss a new mechanism of dust acceleration that acts in a turbulent 
magnetized medium. The magnetohydrodynamic (MHD) turbulence includes 
both fluid motions and magnetic fluctuations. We show that while the fluid 
motions bring about grain motions through the drag, the electromagnetic fluc- 
tuations, can accelerate grains through resonant interactions. In this paper we 
calculate grain acceleration by gyroresonance in Cold Neutral Medium. We con- 
sider both incompressible and compressible MHD modes. We show that fast 
modes dominate grain acceleration. For the parameters chosen, fast modes ren- 
der supersonic velocities to grains that may shatter grains and enable efficient 
absorption of heavy elements. Since the grains are preferentially accelerated with 
large pitch angles, the supersonic grains get aligned with long axes perpendicular 
to the magnetic field. 

Subject headings: dust, extinction-ISM: particle acceleration-kinematics and 
dynamics-magnetic fields 

1. INTRODUCTION 

Dust is an important constituent that is essential for heating and cooling of the inter- 
stellar medium (ISM). It interferes with observations in the optical range, but provides an 
insight to star-formation activity through far-infrared radiation. It also enables molecular 
hydrogen formation and traces the magnetic field via emission and extinction polarimetry 
(see Lazarian 2003). The basic properties of dust (optical, alignment etc.) strongly depend 
on its size distribution. The latter evolves as the result of grain collisions, whose frequency 
and consequences(coagulation, cratering, shattering, and vaporization) depend on grain rel- 
ative velocities (see discussions in Draine 1985, Lazarian & Yan 2002a,b). 
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It is known that all these problems require our understanding of grain motions in turbu- 
lent interstellar medium. Although turbulence has been invoked by a number of authors (see 
Kusaka et al. 1970, Draine 1985, Ossenkopf 1993, Weidenschilling & Ruzmaikina 1994) to 
provide substantial grain relative motions, the turbulence they discussed was not magnetized. 
In a recent paper (Lazarian & Yan 2002a, hereafter LY02) we applied the theory of Alfvenic 
turbulence (Goldreich & Schridhar 1995, henceforth GS95, Cho, Lazarian & Vishniac 2002a 
for a review) to grain acceleration through gaseous drag. 

Here we account for the acceleration that arises from the resonant interaction of charged 
grains with MHD turbulence. To describe the turbulence statistics we use the analytical fits 
to the statistics of Alfvenic modes obtained in Cho, Lazarian & Vishniac (2002b, hereafter 
CLV02) and compressible modes obtained in Cho & Lazarian (2002, hereafter CL02). 

2. Acceleration of Grains by Gyroresonance 

Turbulent acceleration may be viewed as the acceleration by a spectrum of MHD waves 
that can be decomposed into incompressible Alfvenic, and compressible fast and slow modes 
(see CL02). There exists an important analogy between dynamics of charged grains and 
dynamics of comic rays (see Yan & Lazarian 2002, henceforth YL02) and we shall modify 
the existing machinery used for cosmic rays to describe charged grain dynamics 1 . The energy 
exchange involves resonant interactions between the particles and the waves. Specifically, the 
resonance condition is uj — k\\Vfi = nQ, (n = 0, ±1, ±2...), where u is the wave frequency, fey 
is the parallel component of wave vector k along the magnetic field, v is the particle velocity, 
fi is the cosine of the pitch angle relative to the magnetic field, Q = qB / (mc) is the Larmor 
frequency of the particle. The sign of n denote the polarization of the wave. + represents 
left hand polarization and — is for right hand polarizarion. Basically there are two main 
types of resonant interactions: gyroresonance acceleration and transit acceleration. Transit 
acceleration (n = 0) requires longitudinal motions and only operates with compressible 
modes. It happens when k\\v/i = uj. As the phase speed is Vf ast > Va for fast waves, where 
Va is the Alfven speed, it is clear that it can be only applicable to super- Alfvenic grains that 
we do not deal with here. 

Gyroresonance occurs when the Doppler-shifted frequency of the wave in the grain's 
guiding center rest frame uj gc = uj — k\\Vji is a multiple of the particle gyrofrequency, and the 
rotating direction of wave electric vector is the same as the direction for Larmor gyration 



1 In what follows we assume that the time scale for grain charging is much shorter than the grain Larmor 
period. 
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of the grain. The gyroresonance scatters and accelerates the particles. The efficiency of the 
two processes for charged grains can be described by the Fokker-Planck coefficients and 
Dpp/p 2 , where p is the particle momentum (see Schlickeiser & Achatz 1993, YL02). The 
ratio of the two rates depends on the ratio of the particle velocity, the Alfven speed and 
pitch angle, p 2 D fJifl / D pp = [(v(/Va) + At] 2 , where £ = 1 for Alfven waves and ( = k\\/k 
for fast modes. We see that the scattering is less efficient for sub-Alfvenic grains unless 
most particles move parallel to the magnetic field. We shall show later that as the result of 
acceleration, p will tend to 0. Therefore in the zeroth order approximation, we ignore the 
effect of scattering and assume that the pitch angle cosine p does not change while being 
accelerated. In this case, the Fokker-Planck equation, which describes the diffusion of grains 
in the momentum space, can be simplified (see Pryadko & Petrosian 1997): 



df 9/ 1 9 2n , ,9/ 

m +v ^ = ^d-p pDM W {) 

where / is the distribution function. Apart from acceleration, a grain is subjected to gaseous 
friction. For the sake of simplicity we assume that grains are moving in respect to stationary 
gas and turbulence provides nothing but the electromagnetic fluctuations. The acceleration 
of grains arising from the gas motion is considered separately (see LY02). We describe the 
stochastic acceleration by the Brownian motion equation: mdv/dt = —v/S + Y, where m is 
the grain mass, Y is the stochastic acceleration force, S = tdrag I'm is the mobility coefficient, 
tdrag is the gas drag time (Draine 1985, LY02). The drag time due to collisions with atoms 
is essentially the time for collisions with the amount of gas with the mass equal to that of 
the grain, t drag = (ap gr /n n )(7r/8m n kBT) 1 ^ 2 , where a is the grain size, T is the temperature, 
Pgr is the mass density of grain. We adopt p gr = 2.6gcm -3 for silicate grains. The ion-grain 
cross-section due to long-range Coulomb force is larger than the atom-grain cross-section. 
Therefore, in the presence of collisions with ions, the effective drag time decreases by the 
factor given in Draine & Salpeter (1979). When the grain velocity gets supersonic (Purcell 
1969), the gas drag time is given by t s drag = t drag /(0.75 + 0.75c 2 /v 2 - c 3 s /2v 3 + c 4 s /5v A ), where 
c s is the sound speed. 

Multiply the Brownian equation above by v and take the ensemble average, we obtain 

d < v 2 > < v 2 > 

m ^r- = + < £ >• (2) 

Following the approach similar to that in Melrose (1980), we can get from Eq.(l) the energy 
gain rate < e > for the grain with pitch angle p 
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1 d 

< e >= -tHVA^))- (3) 

where the Fokker-Planck coefficient D pp (fi) is given below. 

Adopting the result from quasi-linear theory (hereafter QLT, see Schlickeiser & Achatz 
1993), the momentum diffusion coefficient is (see YL02) 



vpy ' 2w 2 J T ,; 2 + (uj-k\\vii-nn) 2 

+ ^-i(^)^(k)) 

" Jra+l( lfiT )Jn - l( lfir } 

(e^ £ (k) + e- i2 ^(k))], (4) 

where is the nonlinear decorrelation time and essentially the cascading time of the tur- 
bulence, k± = ^k 2 x + ky is the perpendicular component of the wave vector, v± is the per- 
pendicular component of grain velocity, = tan~ 1 (k y /k x ), K a p(k) is the velocity correlation 
tensor and will be given in §3, C,1Z = (x ±iy)/^/2 represent left and right hand polarization. 
At the magnetostatic limit, — > oo, the so-called Breit-Wigner-type function transfers into 
5 function, i.e., t^ 1 /{t^ 2 + (uj — k\\Vfx — nQ) 2 ) — > tx5(uj — k\\vfj, — nf2). Magnetostatic limit 
is correct for fast moving particles (see YL02), but for sub-Alfvenic grains we should use 
Eq.(4). However, we should not integrate over the whole range of fey, because the contri- 
bution from large scale is spurious (see YL02). This contribution stems from the fact that 
in QLT, an unperturbed particle orbit is assumed, which results in non-conservation of the 
adiabatic invariant £ = mv±/2B , where B is the large-scale magnetic field. Noticing that 
the adiabatic invariant is conserved only when the electromagnetic field varies on a time scale 
larger than we truncate our integral range, namely, integrate from k res instead of the 

injection scale L~ l . For Alfvenic turbulence uj = |A;|||Va, the resonant scale corresponds to 
|&||,res| — \£1/(Va — v/jlC) |. For fast modes in a low (3 medium (where (3 = P ga s/Pmag = 2c 2 s /Va 
is the ratio of gaseous pressure and magnetic pressure), u = kVA, the resonant scale is 
k res = {Va — v[/,£)\. The upper limit of the integral k c is set by the dissipation of the 
MHD turbulence, which varies with the medium. 

Integrating from k res to k c , we can obtain from Eq.(4) and (3) the energy gain rate 
< e > as a function of v and /i. Then with the < e > known, we can estimate the grain 
acceleration. Solving the Eq.(2) iteratively, we can get the grain velocity as a function of 
time. We check that the grain velocities converge to a constant value after the drag time. 
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Thus inserting the acceleration rate by fast and Alfven modes into Eq.(2), we can obtain 
the final grain velocities as a function of //. As < e > increases with pitch angle, grains gain 
the maximum velocities perpendicular to the magnetic field and therefore the averaged \x 
decreases. 



3. MHD Turbulence and Its Tensor Description 

Unlike hydrodynamic turbulence, Alfvenic turbulence is anisotropic, with eddies elon- 
gated along the magnetic field. The Alfvenic turbulence is described by GS95 model which 
postulates that k±Vk ~ k\\V A - This may be viewed as coupling of eddies perpendicular to 
the magnetic field and wave-like motions parallel to the magnetic field. For magnetically 
dominated, the so-called low (3 plasma, CL02 showed that the coupling of Alfven and com- 
pressible modes is weak and that the Alfven and slow modes follow the GS95 spectrum. This 
is consistent with the analysis of HI velocity statistics (Lazarian & Pogosyan 2000, Stan- 
imirovic & Lazarian 2001) . According to CL02, fast modes are isotropic. In what follows, 
we consider both Alfven modes and compressible modes in low f3 plasma. 

Within the random-phase approximation, the velocity correlation tensor in Fourier space 
is (see Schlickeiser & Achatz 1993) < v a (k, t)v* p (k', t + r) > jV\ = 5(k - k')K a/3 (k)e- T / Tk , 
where v a ^ is the time-dependent velocity fluctuation in k space associated with the turbu- 
lence. 

The velocity correlation tensor for Alfvenic turbulence is (CLV02), 



T -1/3 

r k = (L/V A )(k ± L)~ 2 ^ ~ (A^r 1 (5) 

where Iy = {Sij — kikj/k 2 } is a 2D tensor in the x — y plane which is perpendicular to the 
magnetic field, L is the injection scale, V is the velocity at the injection scale. Velocity 
fluctuations related to slow modes are subdominant for magnetically dominated plasmas 
(CL02) and we do not consider them. 

Fast modes are isotropic and have one dimensional energy spectrum E(k) oc k~ 3 ^ 2 
(CL02). In low P medium, the velocity fluctuations are always perpendicular to B for all 
k, the corresponding correlation is (YL02) 

r -i/2 

A«(k) = — — Jijk~ 7/2 , r k = (k/L)' 1 / 2 x V A /V 2 (6) 
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where Jij = kikj/k\ is also a 2D tensor in x — y plane. 

4. Results 

We consider a typical cold neutral medium (CNM), T = 100K, n n = 30cm -3 , B = 
6.3/zG. Here, we only consider large grains (lCT 6 cm< a < lCT 4 cm), which carry most grains 
mass (~ 80%) in ISM. The mean grain charge was obtained from the average electrostatic 
potential < U > in Weingartner & Draine (2001). MHD turbulence requires that fluid 
velocities are smaller then the Alfven speed. Therefore we assume that the injection of energy 
happens at the scale L where the equipartition between magnetic and kinetic energies, i.e., 
V = Va, is reached. We assume that the velocity dispersion at the scale / = lOpc is 5km/s 
and that the turbulence at large scales proceeds in tenuous warm media with Alfven speed 
larger or equal to 5km/s. In partially ionized medium, a viscosity caused by neutrals results 
in decoupling on the characteristic time scale (see LY02) 

tdamp ~ ^n 1 ^" 2 ~ (^n)"^" 2 , (7) 

where v n is the kinetic viscosity, l n is the neutral mean free path, v n is the thermal velocity 
of neutrals. Given the parameters above, /„ = 7 x 10 12 cm. When its cascading rate 
t^ 1 = h\\VA ~ k 2 / 3 L^^Va (see Eq.5) equals to the damping rate 2 t~£ mp (k±), Alfvenic 
turbulence is assumed damped 3 . This defines the cutoff wave number of the turbulence 
fc|| jC = 2.4 x 10~ 16 cm~ 1 and the time scale r c = (/c|| iC Va) _1 = 1.7 x 10 10 s. Assuming that 
the grain velocities are smaller than Alfven speed, we can find that the prerequisite for the 
gyroresonance |^|| 5 c| ^ |^||res| Is the same as r L = 2tt/\Q\ > r c , the condition for effective 
hydro drag (see LY02). Thus we see that Alfven modes cannot accelerate grains (with 
a < 2 x 10~ 5 cm) through gyroresonance unless the velocities of these grains are already 
super-Alfvenic. The cutoff of fast modes corresponds to the scale where the cascading time 
scale Tfc = tdamp and this gives the cutoff wave number k c = 4.9 x 10 _15 cm _1 . In the present 
paper we consider neutral gas of low ionization, and therefore the damping owing to ions, 
including collisionless damping, is disregarded (compare to YL02). Using the procedure 
described in §2, we obtain the grain velocities for different pitch angles. The acceleration is 
maximal in the direction perpendicular to the magnetic field. Those values are presented in 



2 As pointed out in §2, the eddy motions happen in the direction perpendicular to the magnetic field, so 
fc_i_is used to calculate tdamp for Alfven modes. 

3 Thus we ignore the effect of slowly evolving magnetic structures associated with a recent reported new 
regime of turbulence below the viscous damping cutoff (Cho, Lazarian & Vishniac 2002c). 
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Fig.l. If averaged over //, the velocities are smaller by less than 20%. 

In order to compare different processes, we account for the hydro drag (see LY02 for a 
discussion of Alfven-induced drag). The fast modes cause the relative movement of the grain 
to the ambient gas also by gaseous drag. Unlike gyroresonance, this relative motion arises 
from the decoupling from the gas. At large scale grains are coupled with the ambient gas, 
and the slowing fluctuating gas motions will only cause an overall advection of the grains 
with the gas (Draine 1985), which we are not interested. The largest velocity difference 
occurs on the largest scale where grains are still decoupled. While in the hydrodynamic case 
the decoupling happens on the time scale tdrag, the grains are constrained by the Larmor 
gyration in MHD case unless td rag < t~l (see LY02). The latter condition is true for high 
density gas, however (see Lazarian & Yan 2002b). The velocity fluctuations for fast modes 
scales as v k oc /c~ 1//4 oc cj -1 / 4 , where u = Wa is the frequency of fast modes (see Eq.(6)). In 
low (5 medium, the velocity fluctuations are perpendicular to the magnetic field. Therefore, 
grains get velocity dispersions perpendicular to the magnetic field, v ~ V(7x/r. ma , E ) 1 / 4 , where 
T~max = L/V is the time scale at the injection scale of turbulence. Then we also need to 
consider the effect of damping. Similar to Alfven modes, the condition for effective hydro 
drag tl > r c = 2it/uj c = 5 x 10 9 s is the same as k c > k res for sub-Alfvenic particles, 
which is the requirement for gyroresonance. Our calculation shows that the corresponding 
grain size is 4 x 10 _6 cm. For smaller grains, their velocities are reduced v ~ v c x tl/t c ~ 
U(r i /r m , aa; ) 1 / 4 (rL/r c ) 3 / 4 , where v c is the velocity of turbulence at the damping scale (see 
Lazarian & Yan 2002b). In Fig.l we plot the velocity of grain with \x = as a function of 
grain size since all the mechanisms preferentially accelerate grains in this direction. 

How would the results vary as the parameters of the partially ionized medium vary? 
From Eqs.(2), (3)&(4), we find the grain velocity is approximately equal to v ~ (< e > 
tdra 9 /m) 1/2 ~ 2.5 x 10 4 L^ /4 V 5 Bl, /2 (q e / al)°H° d f ag cm/s, where L 10 = L/lOpc is the injection 
scale defined above, V 5 = V^/Skm s _1 , = B/lfxG, q e = g/lelectron. Noticing that the 
hydro drag by fast modes decreases with the magnetic field, v oc B~ 1 ^ beyond the damping 
cutoff and v oc B^ 1 below the cutoff (see the last paragraph), we see that the relative 
importance of the gyroresonance and hydro drag depends on the magnitude of the magnetic 
field. 

It has been shown that the composition of the galactic cosmic ray seems to be better 
correlated with volatility of elements (Ellison, Drury & Meyer 1997). The more refractory 
elements are systematically overabundant relative to the more volatile ones. This suggests 
that the material locked in grains must be accelerated more efficiently than gas-phase ions 
(Epstein 1980; Ellison, Drury & Meyer 1997). The stochastic acceleration of grains, in this 
case, can act as a preacceleration mechanism. 
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Fig. 1. — Grain velocity vs. size owing to different acceleration processes in CNM. The 
solid line represents the gyroresonance with fast modes. The dashdot line refers to the 
gyroresonance with Alfven modes. Gyroresonance with both modes works only for large 
grains owing to the cutoff by viscous damping. The cutoff scales for fast and Alfven modes 
are different mostly due to the anisotropy of Alfven modes. The dotted line is the result 
from hydro drag with Alfven modes (see LY02), the dashed line represents the hydro drag 
with fast modes. The corresponding grain Larmor frequency Q and are also given along 
the abscissas. 
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Grains moving supersonically can also efficiently vacuum-clean heavy elements as sug- 
gested by observations (Wakker & Mathis 2000). Grains can also be aligned if the grains get 
supersonic (see review by Lazarian 2003). Indeed, the scattering is not efficient for slowly 
moving grains so that we may ignore the effect of scattering to the angular distribution of 
the grains. Since the acceleration of grains increases with the pitch angle of the grain (see 
Eqs.(3) and (4)), the supersonic grain motions will result in grain alignment with long axes 
perpendicular to the magnetic field. 

It is believed that silicate grains won't be shattered unless their velocities reach 2.7km/s 
(Jones et al. 1996). Nevertheless, the threshold velocities depend on the velocity of turbu- 
lence at the injection scale and on the grain structure, i.e., solid or fluffy. Thus shattering 
of the largest grains is possible. 

5. SUMMARY 

In the paper we showed that 

1. Fast modes provide the dominant contribution for the acceleration of charged grains. 
The velocities obtained are sufficiently high to be important for shattering large grains and 
efficiently absorbing heavy elements from gas. Alfven modes are not important because of 
their anisotropy. 

2. Depending on the relative importance of the magnetic field, gyroresonance (strong 
B) or hydro drag (weak B) by fast modes dominates grain acceleration. For small grains, 
hydro drag by fast modes is the most important while the gyroresonance is not present as 
the turbulence at the resonant frequencies gets viscously damped. 

3. In low (3 medium, all the mechanisms tend to preferentially accelerate grains in the 
direction perpendicular to the magnetic field. Among them, gyroresonance with fast modes 
can render grains with supersonic motions, which can result in grain alignment perpendicular 
to the magnetic field. 

We acknowledge stimulating communications with E. Zweibel, which were very impor- 
tant at the initial stage of this work. We thank J. Mathis, J. Cho and R. Schlickeiser for 
helpful discussions. This work is supported by NSF grant AST0 1-25544. 
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